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ABSTRACT 



The sensitivity of planetary waves to various initial 
conditions and thermal forcing is examined using a linear 
global primitive equation spectral model. 

Initial conditions are obtained by switching on a 
analytic heat source and then integrating the model equa- 
tions out to 30 days. An averaging procedure is used to 
eliminate any transient modes which remained after integra- 
tion so that the initial conditions represent steady state 
solutions. Additional integrations are performed in which 
'errors’ are introduced into the forcing and initial 
conditions. 

Results of the study are examined using polar phase vs 
amplitude plots (harmonic dials) of various spherical 
harmonics. Results indicate that planetary waves are not 
sensitive to errors in the forcing or initial conditions. 
However, this laclc of sensitivity is most likely due to the 
simplified initial conditions used in the model rather than 
to the inherent sensitivity of planetary waves. 
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I. INTRODUCTION 



Atmospheric predictability studies (Lorenz, 1969) indi- 
cate each scale of motion has a limit to it’s predict- 
ability. Small scale motions are theoretically predictable 
to an hour, synoptic scales to a few days and planetary 
scales to a few weeks. However a number of studies, Lambert 
and Merilees (1978), Baumhefner and Downey (1978), and Morse 
(1983), have shown that synoptic scale motions, rather than 
planetary scale motions, are the most accurately forecast. 
Fig. 1.1 is taken from Daley et al (1981). It gives the 
500mb geopotential forecast error from the 6-level National 
Center for Atmospheric Research global circulation model as 
a function of wavenumber- The error for each wavenumber has 
been normalized by the natural atmospheric variance for that 
wavenumber. Note that as previously stated the synoptic 
scales (wavenumbers 4-8) are in fact the most accurately 
forecast. This result does not seem to be dependent on the 
type cf model used. Baumhefner and Downey examined a number 
of different models and found similar results- The error in 
forecasting planetary scale waves might not seem to be that 
important given that the error in the synoptic scales is 
smaller and that these are the scales of motion which 
produce most of the day to day weather changes- The impor- 
tance of an accurate planetary scale comes to light when one 
starts to consider medium to long range forecasts (up to 1 0 
days). Since the planetary waves often act to steer the 
smaller synoptic scale disturbances an improved planetary 
scale forecast would presumably lead to an improved forecast 
on the synoptic scale. In addition, the planetary waves 
contain a major portion of the eddy kinetic energy in the 
atmosphere and again a better long wave forecast would lead 
to a better synoptic scale forecast. 
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There are at least two possible reasons why planetary 
waves are not forecast as well as theory suggests. One is 




Figure 1. 1 Forecast Errors for 1, 2 
and 3 Days of Integration. 

that the model dynamics for planetary waves are inadeguate. 
Given that planetary waves are guasi-stationary in nature 
and are forced in part by differential heating and orography 
it would seem reasonable to suspect that inadequacies exist 
in the dynamics of the model's forcing. Another possibility 
is that current initialization procedures are not adequate 
for planetary scales. 

It is the purpose of this study to examine the nature of 
planetary waves and to deternine the sensitivity of 
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planetary waves to icitial conditions and forcing. Chapter 
two of this thesis describes the model used in this study. 
The third chapter outlines the nature of the experiment. In 
the fourth chapter the results of the experiment are 
presented and in the final section these results 
discussed and some conclusions are drawn. 
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II. MODEL DESCRIPTION 



The model used in this study is a baroclinic spectral 
transform model which was developed by Rosmond (1977). 
Similar formulations have been carried out by Hoskins and 
Simmons (1975) and by Bourke ( 1974). The model is config- 
ured to include friction and/or diabatic heating- The 
specifics of how friction and diabatic heating are included 
in the model will be discussed in the next section. The 
basic equations of the model, in sigma coordinates, are as 
follows: 



1^ = -V*(C + f)^-)<-^x(RT^q + 6|^)+ic*VxF 

d L p G 

Iy = k*Vx( ; + f )V-V- (RTVq+o|^)-V^ (4)+j ) +V • F 



^ - 
8t ■ 



-D-V* Vq- 



da 

da 



ii 

dt 



_V*V9+d|^+ 

o O 



Q_ 

PCp 



da ~ a 



where : 

^ - vorticity 
D - divergence 
T - temperature 
0 - potential temperature 
TT - surface pressure 
V - horizontal velocity vector 
d - geopotential height 
R - gas constant 

Cp - specific heat at constant pressure 
f - Coriolis parameter 



( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 
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o - vertical coordinate (o = P/tt) 

0 - vertical velocity (d = ^ ) 
q - In 

P - P 
< - R/Cp 

F - frictional force 

The continuity equation (Eq. 2.3) may be rewritten by 
inteqrating with respect to signa and applying the boundary 
conditions d(0)=d{1)=0. Thus the integral of Eg. 2.3 maybe 
written as 



fa = .d*g <2.6) 

where ( ) = /^( )do and G= V *Vg. 

The vertical velocity, d/ obtained diagnostically 

by substituting Eg. 2.6 into Eg. 2.3 and integrating in the 
vertical to obtain 



d = (D+G)a-/°(D+G)do ^ ^ ^ 

which uses d (O) =0 

The first law of ther mod yna nics Eg. 2.4 can be written 



at 



-V-VT-o‘^d|^(Ta"'^)+KT(|^+V-^q)+^ 



( 2 . 8 ) 



In order to apply semi-implicit differencing it is necessary 
to divide temperature as follows: 



T = T*(o)+T'(o, X, T, t) 

where T* represents an appropriately averaged temperature. 

The basic equations can be conveniently written in 
spherical coordinates by defining the following operator; 
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( 2 . 10 ) 



ct(a,b) 



1 3a 3b 
3X 3x 



Using Eqs. 2.9 and 2.10 the basic equations can be written 
as follows: 



II: = -a(A,B) 

d L 



= -a(A,B)4^(E+4)+RT*q) 



3t 



-a(UT' ,VT’ )+r7T'-a^a|-(Ta ’^)+kT(G-G-D) 
’ 3a 



3t 



^ = -(D+G) 






where 



„ , .3U xRTl3q „ 



Q z y_ is. +v^ 

l-x2 3X 3x 



E = (U^+V^)/2(l-x^) 



U = u(cos'l')/r 
V = v(cos'y)/r 
X = sin'}' 
f = latitude 



(2. 11) 
( 2 . 12 ) 
{2, 13) 
(2. 14) 
(2. 15) 



X = longitude 



1 4 



Eys. 2.10 - 2.15 are the basic equations used in the model. 
These equations are represented spectrally in the horizontal 
and finite differenced in the vertical. 

A. VERTICAL STRUCTOEE 

The vertical structure of the model follows the develop- 
ment given by Arakawa and Suare 2 (1983). The variables are 
staggered in a so that K, D, U , V and T are carried at the 
mid-point of each layer where ^ ^ is carried at the 
top and bottom of each layer where a . The vertical 
structure is illustrated in Fig. 2.1. The finite difference 
form of Egs. 2.10 - 2.15 are 




(2. 17) 



(2. 16) 



(2. 18) 




(2. 19) 




( 2 . 20 ) 



K 




( 2 . 21 ) 




k 

(G+D)-^ (G.+D.)Ao. 

3=1 3 3 3 



( 2 . 22 ) 



where 
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^,= 

D, T, Aa, 0, 









# 



^ K -1 ^ K -1 *^<-1 ^ K -1 ^ K -1 ^ K -1 

t D 1 Ag 6 

5 K K K K r K 






(7 



K + 1 





I 




Figure 2. 1 Yertical Structure. 
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-(CpRTj;./r)|^ +(/r^Vr)F^ 
G = Vj^*^7r/ir 



D = ^*V. 



k 





Egs. 2.16 - 2.22 can be written in matrix form, such 

that the terms on the right hand sides contain all the terras 
which are to be evaluated explicitly and the left hand side 
contains those terms which are to be evaluated implicitly. 
Eg. 2.20 can be combined with an integrated finite differ- 
ence form of Eg. 2. 15 to obtain 

(j) = C 

where ^ is a sguare matrix and the other quantities are 
column vectors. The finite difference form of the surface 
pressure tendency eguation(Eg. 2.14) is 




(2. 24) 



which can be written in matrix form as 



|| = -n'^(G + D) 

O L 



(2.25) 
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T . 

where N is the transpose of a constant column vector. 
Similarly E^. 2.22 can be written 



o = Z(G+D) 

Ihe next to last term in Eg. 2- 19 is 

k\- 1 < : 



(2. 26) 



) 



(2.27) 



For the purpose of se ni-implicit formulation, the tempera- 
ture is separated according to Eg. 2.9. The mean part of 
that term can be written as 



= M(G+D) 

K+ 1 )C_ 1 

Egs. 2.17, 2.19, 2.23, and 2.25 may now be written: 



(2. 28) 



l^+V^ ($ ' +RT*q)=a(B,-A)-V^E = 


(2. 29) 


I^+n'^D = -N^G 

O L 


(2. 30) 




(2.31) 


II 

In 


(2. 32) 



where <T*N and <})= ‘i' ~ 'f’g and and K represent 

terms which have not been explicitly separated out. 

The semi-implicit time differencing is achieved by eval- 
uating the terms on the left hard sides of Egs. 2.29, 2.30 

and 2.31 implicitly. The remaining terms and Eg. 32 are 
evaluated explicitly using leapfrog differencing. The 
difference eguations can now be written 

D .,+AtV^(C T = 

^ = n + 1 — n+ 1 



n + 1 



D -,-AtV^(C^ ^) + 2At(K,^) 

= n-1 _ ^n-1 



{2. 33) 



n- 1 



D n 
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(2. 34) 



q„,j.AtSTB = 

Now the following equation for D can be found by substi- 
tuting Eg- 2.34 and 2.35 into 2.33: 

-RT'-N G 
= n 

where the matrix operator B is 



= At^ (C Q±RTn'^)V^+^ 

B. SPECTBAL FOBMDLATION 



The equations (2. 10} - (2. 1 5) are represented spectrally 
in the horizontal. The variables are represented as 
follows: 



C(X. X, o, t) = i___j t)p”(x)e“" 

J J 

V V ^ni^rn 

L T b 1 I G Y 

m=-J n= m n n 



(2. 38) 



m 



m -m 



where C is some variable and (Cn ) = (~1) . Cn 

m = zonal wavenumber 

n = meridional index, and n-jm) gives the number of zeros 
between the poles (-1 < x < 1) of the associated Legendre 
function 

J = truncation limit (for this study triangular truncation 
was used with J=3) 



(1-1) /2 nondimensional zcnal coordinate index 
(1 < 1 < 16) 
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Note that the separation is such that the 
are functions of time and the vertical and 
harmonic Y^are horizontal functions of space 
zation and orthogonality properties of the 
coefficients to be obtained as follows: 



coefficients 

the spherical 

. The normali- 

Y’’^ allow the 
n 



"n 



<c,y"'> 

’ n 



1 , 2i\ 

4ti 



/ t^CY^^dxdX 
-1 n 



(2.39) 



The non-linear terms are computed using the transform method 
following Kaltiner and Williams (1980). The longitudinal 
direction is treated with a Fast Fourier Transform and the 
latitudinal direction uses Gaussian Quadrature. The number 
of latitudes, N and longitudes H satisfy. 

N > 3J/2 +1 M > 3J+ 1 



The number of points are chos 
ciliasing from the product term 
H=16. It was discovered j 
completed that number of Gaussi 
the latitudinal integrals was a 
However, model integrations u 
qualitatively different than th 
Even the quantitative differenc 
mostly in the non-linear inte 
examined in this study are car 
of the model, one would not exp 



en so that there will be no 
s. For this study N = 4 and 
ust before this study was 
an latitudes used to compute 
ctually J+1 and not 3J/2 + 1. 
sing J +1 latitudes are not 



cse using 


3J/2 + 1 


latitudes. 


€s are small, and 


will occur 


grations. 


Since 


the cases 


ried out 


using a 1 


inear form 


€ct much 


error. 
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III. EX£MIMMT design 



The intent of this study is to specify a known analytic 
heating function on a planetary scale and then integrate the 
equations until the model atmcsphere reaches a steady or 
quasi-steady state. These steady state solutions will be 
altered and/or the heating will be changed. The equations 
will then be integrated again to determine the sensitivity 
of planetary waves to changes in the forcing and initial 
conditions 

Heating is introduced into the model via the thermody- 
namic equation (2.4). Q is specified to be a Legendre 
function of the form: 



Q(X, X, o) = i z I I A"'(a)p"'(x)e^^^ 
m=-J n= m n n 



(3.1) 



is specified to have the vertical structure: 



n 




(3.2) 



Only one wave form is forced at a time and the other heating 
amplitudes are set to zero. The constant A is chosen to 
give a reasonable temperature response. In most cases a 
temperature response of a few degrees is used. 

Before the heating is '*turEed on" the model atmosphere 
is specified to be in a state of solid rotation 

^ = cos 
V = 0 

such that U = 20 m/s(15.5 °/day) at the equator. The 
surface pressure field is chosen to balance the mean flow. 
The initial vertical temperature distribution is given by 
the standard atmosphere values consistent with the pressure 
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distribution given by solid rotation. A lapse rate of 6.5° 
per kilometer is specified to the top of the model 

atmosphere. 

TABLE I 

Haves Forced for Linear Cases 

Case Have F or ced 

m n 

1 1 1 

2 1 2 

3 1 3 

4 2 2 

5 2 3 

6 3 3 

In all of the experiments the model is run in a linear 
mode. In this mode the amplitudes of all zonal wavenumbers 
except those being forced are held constant. The model is 
run in a linear mode in these experiments to allow a more 
definitive interpretation of the results. Six linear cases 
are integrated, each case with a different forced wave 
form. For example in case 1 (see Table I) zonal wavenumber 
one(m=1), meridional index one (n = 1) is forced. The results 
from these cases are fairly similar so only the longest wave 
case (case 1) is analyzed in detail. 

Two classes of cases are examined using the steady state 
values. In the first, the magnitude of the heating function 
is changed and the model is integrated from the steady state 
values. In the second, the nagnitude of the steady state 
values are changed and then used as initial conditions for 
further model integrations. In the these cases the heating 
function is the same as the one originally used to bring the 
model to a steady state. 
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17. BESDLTS 

Due to the similarity of the six linear cases only cases 
1, 5 and 6 will be discussed. These three cases represent 
the three zonal wavenumbers ( 1 , 2 and 3) and the two basic 
meridional structures possible in the model. Cases 1 and 6 
are symmetric in their meridional structure while case 5 is 
anti-symmetric. The primary tool used to analyze these 
cases was the harmonic dial. These diagrams are polar plots 
of phase and amplitude for a given spherical harmonic. This 
plot has the advantage that the guasi-stationary and tran- 
sient nature of the waves are easily seen. Fig. h. 1 is a 
good example upon which the merits of the harmonic dial can 
be explained. In this figure the phase and amplitude for 
the temperature wave for case 1 are plotted every 12 h. The 
elapsed time since the start of the integration is indicated 
every 5 days. For an ideal case where there is a stationary 
and a transient component one would expect to see a circular 
pattern which is confined to a particular quadrant. This is 
because as the transient wave moves through the stationary 
wave it will come in and out of phase with the stationary 
wave. The maximum amplitude occurs when the transient and 
stationary waves are completely in phase and the minimum 
amplitude occurs when the waves are 180° out of phase. A 
rough estimate of the phase speed of the transient wave can 
be obtained by observing the change in phase with time as 
plotted on the diagrams. The direction of propagation can 
be determined by observing the direction of the phase change 
on the dial. A counterclockwise change in phase indicates 
eastward propagation, a clockwise change indicates westward 
propagation. In Fig. 4. 1 the phase speed of the temperature 
wave at a = .925 can be easily estimated. The period of the 
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Have is 21.5 days# 
is 360/21.5=16.74. 



thus the phase speed in degrees per day 
The direction of propagation is east- 



PHASE VS AMPUTUDE 
FOR TEMPERATURE 




Figure 4.1 Harmonic Dial for Case 1 Temperature Component. 

ward. Another important point to note here is that, given 
the uniformity of the change in phase, the phase speed of 
the transient component is constant. The harmonic dials for 
surface pressure, temperature and the zonal(U) and meri- 
dional (V) wind components for the lowest model level are 
shown in Figs. 4.2 - 4.7 . Only those components which have 
significant amplitude are shown in these figures. In all 
cases the pressure and winds waves show a much more irreg- 
ular pattern than the temperature wave. These irregulari- 
ties might be due to the presence of a number of transient 
modes components with periods much less than that indicated 
by the overall pattern. Table II gives the estimated period 
(in days) and phase speed (in °/^ays) for the temperature, 
surface pressure and the V and U wind component waves for 
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PHASE VS AMPUTUDE 



PHASE VS AMPUTUDE 
FOR SURFACE PRESSURE FOR SURFACE PRESSURE 





PHASE VS AMPUTUDE 
FOR TEMPERATURE 




J 



Figure 4.2 Harmonic Dials for Temperature 
and Surface Pressure for Case 1. 
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PHASE VS AMPLITUDE 



PHASE VS AMPUTUDE 




Figure 4.3 Harmonic Dials for Meridional (0) 
and Zonal (D) Wind Case 1. 
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PHASE VS AMPUTUDE 



FDR SURFACE PRESSURE 




PHASE VS AMPUTUDE 
FOR TEMPERATURE 




Figure 4.4 Harmoaic' Dials for Temperature 
and Surface Pressure for Case 5. 



27 



PHASE VS AMPLITUDE 



FOR V COMPONENT OF WIND 




PHASE VS AMPLITUDE 



PHASE VS AMPLITUDE 



FOR U COMPONENT OF WIND 




FOR U COMPONENT OF WIND 




Figure 4.5 Harmonic Dials for Her idional (0) 
and Zonal (0) iind Case 5. 
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PHASE VS AMPLITUDE 




PHASE VS AMPLITUDE 
FOR SURFACE PRESSURE 




Figure 4.6 Same as Fig 4.4 for Case 6. 
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PHASE VS AMPLTUDE 



FOR V COMPONENT OF WIND 




PHASE VS AMPLfTUDE 
FOR U COMPONENT OF WIND 




Figure 4.7 Same as Fig 4.5 for Case 6 



the three cases being considered. These estimates were 
obtained from Figs. 4.1 - 4.7 . Phase speeds were not 

calculated for the pressure and wind waves for case 5 due to 
the highly irregular nature of those dials. The phase speed 
of the temperature, pressure and wind waves are the same for 
wave component (3,3) . Such is not the case for the other 
cases, though the phase speed are similar. The computed 
phase speed of the different variables for case 1 varies 
from 15.6 ®/day for P to 18.5 ®/day for U. The other cases 
have similar phase speeds for the different components. The 
general pattern which emerges from the data is that the 
phase speed of the transient components seems to be indepen- 
dent of wavenumber and that the waves are moving at approxi- 

TABLE II 

Period (p) , and Phase Speed (c) for Transient Components 



Case Case Case 

“T“ ”5" 





p 


c 


F 


C 


p 


c 


p 


23: 5 


15. 6 






7.5 


16.7 


T 


21.5 


16.7 


12.5 


16.7 


7.5 


16.7 


U 


19.5 


18.4 






7.5 


16.7 


V 


21. 5 


16.7 




❖ 


7.5 


16. 7 



mately the speed of the mean wind (15.5 °/day) . To test this 
hypothesis additional integraticns were accomplished holding 
all other factors the same but increasing the speed of the 
mean wind. The results of these integrations showed an 
increase in phase speed of the component waves corresponding 
to the increase in the mean wind speed 

A. STEADY STATE CASE 

Since the purpose of this study was to see how sensitive 
the steady state solutions are to errors in the forcing and 
initial conditions it is important that the steady state 
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solutions of the model be examined to insure that they are 
indeed steady and that these solutions are consistent with 
known theory. 

To obtain steady state solutions the basic equations 
were integrated up to 30 days. It is obvious from Figs. 
4. 1-4. 7 that after a rather lengthy integration the model 
atmosphere has yet not reached a steady state and that in 
fact the transient modes identified in the last section are 
clearly evident throughout the period of integration. 
Though there are indications that the solutions are slowly 
converging as indicated by the begining of spiral patterns 
in Figs. 4. 1-4. 7 a final steady state would not be reached 
for some time. However, by averaging over the period of the 
transient wave, it’s effect can be averaged out. The 
remaining values are the steady state values. The above 
procedure was performed on case 1 . To insure that the part 
of the solution which remained after averaging was truly the 
steady state solution another integration was performed 
using the averaged values as the initial conditions. 
Harmonic dials were again constructed for each of the vari- 
able wave components. If the fields obtained by averaging 
were truly the steady state solutions, one would expect to 
see no phase or amplitude change with time, i.e. the 
harmonic dial would collapse to a single point. Figs. 4.8 
and 4.9 contain the harmonic dials from a 20-day integration 
initialized with the averaged values for case 1. The 
temperature wave is almost completely stationary, but the 
pressure and velocity components have small growth and move- 
ment. The fact that the pressure and wind component waves 
are not completely stationary is another indication that 
there are transient modes in these waves which are not 
contained in the temperature wave. Consequently , the aver- 
aging which was performed over the period of the temperature 
wave did not average out all of the transient parts of the 
pressure and wind waves. 
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PHASE VS AMPUTUDE 
FOR SURFACE PRESSURE 




PHASE VS AMPUTUDE 
FDR SURFACE PRESSURE 




PHASE VS AMPLITUDE 



FDR TEMPERATURE 




PHASE VS AMPUTUDE 




Figure 4.8 Same as Fig 4.2 for Steady 
State Initial Conditions. 
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Figure 4.9 Same as Fig 4.3 for Steady 
State Initial Conditions. 
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It appears from the above integrations that solutions 
obtained by averaging over the period of the transient 
component are very close to steady state values. In order 
to better understand these solutions they will be compared 
to the steady state solutions from a simplified analytic 
model. This model will be derived in the next section. 



B, THE AHALYTIC HODEL 



The 

motion 



predictive 
in Z=-ln (P/P) 



system of equations 
coordinates are: 



for 



planetary scale 




(4.1) 






+ V 



^ (f) ^ • — • 

+zr(z)+zr’ 



Q(Z) 



(4.2) 



V = jkxVcj) 
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D = ^*V 



Z = -p/p 



(p = 4>(Z)+ip' (x, y, Z) 



r 
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9Z^9Z 



+ K<j) ) 



^(_g 

T ^Cp 



+ 



•1 9T 
H 9Z 



= RT/g 



r 



t 



U1K(S_ 

T ^Cp 



+ 



1 

H 9Z 



Eg. 4.1 is the continuity equation. Eg. 4.2 is the first law 
of thermodynamics and Eg. 4.3 is the geostrophic relation. 
These equations were obtained by scaling the general haro- 
clinic equations (Haltiner and Williams, 1980) with: 
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L/a-l e <.01 e ~ 10 0 Ra/f ~1 

Where L is the horizonal scale, R is the Rossby Numberand e 

is a rotational Froude number. Assuming a steady state 

atmosphere (^= o ) solid rotation (U= 03 a cos'!' ) and 

9 1 

droping ther’(Z) term, because it is small. Eg. 4.2 becomes: 

(4.4) 

components of the geostrophic wind 
ar e : 

V - 1 ii' ('^•5) 

2af^ cos'J'sinT 9A 

U - - 1 94> ' (4.6) 

2a2J^sin'y 3'F 

Using Egs. 4.5 and 4.6 the divergence in spherical coordi- 
nates can be written: 



“lx = Q<2> 

The meridional and zonal 
in spherical coordinates 



„ _ 2ncosH' _ 6V 
liOsTnY ■ f 

Now Eg. 4.1 can be written: 



9Z ; 1 9(j) 

9Z " a^ITsinW 9A 

Assuming solutions 



= 0 

of the 



form : 



(4.7) 



(4.8) 



<t>' = <J>(Z)e^^^P^('l') 


(4.9) 


Z = W(Z)e^^^p'”('i') 
n 


(4. 10) 


Q = N(Z)e^^^pJ^('l') 


(4. 11) 


Eg. 4.4 and 4.8 become: 




ojim'J'C Z) +rw( Z) = N(Z) 


(4. 12) 
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0 



(4. 13) 



w -w- 



im'^> 
sin^ 4' 



Eg. 4.12 and 4.13 can be combined to obtain a single equa- 
tion for W 



oja^J^sin^'l'(W"-W' )+rW = N 



(4. 14) 



Eg. 4.14 is just a 
ential equation for 



W = CtWt+C.W„+W 

11 2 2 p 



2nd order n cnhomogenous ordinary differ- 
H whose general solution is of the form: 



Hhere F]^ , 1^ are the solutions to the homogeneous part of 

Eg. 4.14 and W„ is the particular solution. Assuming N (Z) = 

2 ^ 

? the solution to Eg. 4.14 is 



,, „ r+Z,„ r_Z,. -a^Z 

W = C^e ^ + 0 - 2 ^ +Ae 



= («. 15) 



wher e 



= -2 ± /l/4-r/cja2J^sin'F 






a = seme positive number 



A = NQ/(coa^sin^'}'(a'’+a^ )+r) 

There are two possible solution cases for the homogeneous 
part of Eg. 4.14. If gs is positive the solution consists 
of an exponentially growing and decaying part. After 
applying the boundary conditions w(0)=0 and finite energy at 
Z=~, i.e. {p|(})|2 ) = e~^|<j)12 the solution becomes 



A/ -ot^Z Z/2 -q^Z, 
W(Z) = A(e -e e ^ ) 

An expression for $ (z) can 
4.16 into Eg. 4. 13 



(4. 16) 

be obtained by substituting Eg. 
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<I>(Z) = -a^Hsin^ 'i'A{e^^(r-l) +e ^ 



( 1 + ) } 



{4. 17) 



If g 2 is' negative the solution will be oscillatory in 
of the form 

W(Z) = + {“* 

where 

M = /r/o)a2Qsin2'i'- 1/4 

After applying the radiation boundary condition at Z= « 
w (Z) (i) (Z) >0 and the lower boundary condition w = 0 at Z = 

4. 18 becomes 

W(Z) = A{(|i-u)e^''^‘^^^^^-i(-l-aMe‘“^^} 

Ihe corresponding expression for (z) is 

$(Z) = a2j^sin2'l'A{(|i-y)e^''^'^^’''^^-i(-l-aMe~“^^} 

Eqs. 4.16,4.17 and Egs. 4.20,4,21 can be substituted 
Eqs. 4.9 and 4.10 to obtain expressions for (|)and Z. 
taking the real part the g2 > 0 solution are 

<P' = a^Qsin^TA{e^^(r-l) +6"°^ 1 + aM 1 sin(mA ) 

n 

Z = A(e ^ )cos ( mX ) 

n 

The g2 < 0 solutions are 

(}) ' = a^f^ sin^ TAC ^ ( J-sin( yZ + mX ) 

+ ycos ( y Z+mX ) + ( 1 + a^ ) e ^sin( mX ) } P'” 

n 



(4. 
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A /V f -a^Z f Z/2 / r,N . 

Z = Ale -(e cosCyZ) )cos(mA) 



(4. 25) 



+ (e^'^^sin(viZ) )sin(mA) 

If > 0 the response of the atmosphere to heating is 

trapped near the surface and the vertical structure is one 
in which the waves have no tilt with height. If g2 < o the 
perturbation due to the heating will radiate to great 
heights and the waves will tilt westward. In both cases tlie 
steady state represents a balance between advection, adia- 
batic warming (cooling) and diatatic heating (cooling) . The 
longitudinal and latitudinal structure, for each case, will 
be determined by the forcing specified. That is, if waven- 
umber one is specified in the forcing term, the steady state 
solution will have a wave one form. In addition, since the 
winds are geostrophic their structure will be determined by 
the structure of ' . Eg. 4.24 can be written 

4)' = -A{ (je^'^^sinCyZ) -e^^^ycos(yZ) )cos(mA) 

+ (|e^^^cos(yZ)-ye^^^sin(yZ)-(l + aMe"°‘^^) 
which can be written in the form 

(G(Z)cos(nA+H(Z) ) } A( ’!') P^( T ) 

n 

This form allows comparison with the phase and amplitude of 
the model Legendre functions. If the model steady state 
solutions really are steady state one would expect them to 
at least gualitatively agree with the analytic solutions 
derived above. The differences between the numerical model 
atmosphere and the analytic atmosphere are: 

1. A nondimension vertical coordinate Z is used in the 
analytic model while sigma coordinates are used in 
numerical model, 

2. Friction is not included in the analytic model, and 



(4. 26) 

sinCmA ) } P^ 
n 
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3. r(Z) is assumed to be constant in the analytic model 
while it is not in the numerical model. 
r(Z) is not constant in the numerical model because r (z) is 
constant. The fact that r(Z) is not constant given r (z) is 
can be shown by examining the eguation for r (Z) 



r(z) = 






4(3 + i il) 



T 'Cp H az' 

After converting to 2 coordinates Eg. 4.28 becomes 



(4. 28) 



T ^cp az^ 



(4. 29) 



Since the lapse rate is a constant 6.5 °/Km in the model the 
second term in Eg. 4.29 is constant but T decreases with z 
so the first term is not constant thus T (Z) will not be 
constant in the model atmosphere. However, since the 
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Fi 



ure 4.10 Heating 
a) Numerical model 



Vertical Structure 
and (b) Analytic. 
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be equal to 1/(-ln.7). Fig. 4.10b is a plot of 
converted to P coordinates. Note that the plots are very 

similar up to 400mb. Figs. 4.11a and 4.11b are plots of 
G(Z) and H (Z) vs pressure using the following estimate for 
the various parameters in Eq. 4.19: w =(20 m/s)/a, r~ (Z) 

H2N2, K = 10km, N2 = 10~*s-i, a=6.37 x 10^ m and ^ = 7.292 x 

10 -Ss-i . Figs. 4.12 and 4.13 contain plots of the steady 
state model solutions for the neridional (V) and zonal (U) 
wind vs sigma for case 1. Since the winds are geostrophic 
the vertical structure U and 1 can be compared with the 
vertical structure of 4>. Note the striking similarities 
between the two sets of figures. Both amplitude plots show 
a slight decrease in the lower layers and an exponential 
increase near the top of the atiosphere. Also, given Figs. 
4.14, 4.15 and 4.16, it is obvious that the model’s steady 
state solutions have the longitudinal wave one form speci- 
fied in the case 1 forcing term. These figures are hori- 
zontal cross-section for vorticity, meridional velocity and 
temperature near the top and bottom of the model atmosphere. 
Contour intervals are indicated at the top of each cross- 
section. The vertical scale at the bottom of the figures is 
used to represent topography, which is not included in this 
model. It is clear that the steady model solutions are 
qualitatively consistent with the analytic solutions and 
that the method used to determine these values was a sound 
one. The question might now be asked why are the model 
solutions oscillatory in nature? Is it not possible that 
they could be of the exponential decaying type? These ques- 
tions can be answered by examining Eg. 4.19 more closely. 
Given that 

r/coa^sin^4' > (^-30) 
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Figure 4.11 Analytic ¥ertical Structure 
(a) Aiplitude and (b) Phase. 
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Figure 4. 12 Model TJertical Amplitude 
and Phase Structure for Case 1. 
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Figure 4. 13 Ilodel Vertical Amplitude 
and Phase Structure for Case 1. 



45 



f AT HOUR 0 




Figure 4.14 Steady State Longitudinal 
Vorticity Structure for Case 1. 



46 



■T AT HOUR 0 







47 
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Figure 4. 16 
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the solution will be oscillatory. In the model, co is speci- 
fied so that the zonal wind is equal to 20 m/s at the 
equator so that U= w a cos = U=wa at =0 thus co = U/a. Given 
r~H2 H2, H-10 km, N2 =10-»s-i, a = 6.37 x 10& m, Q =7.292 x 
10~ss-i one can try to solve fcr a latitude where the solu- 
tion will no longer be oscillatory. After using the above 
values Eq. 4.30 becomes 



sin¥ = 20.84 (^-31) 

which has no solution. Thus there is no point where the 
solution would not be oscillatory. To obtain non- 

oscillatory solutions one would have to increase the zonal 
wind speed to almost 100 m/s. Using these steady state 
solutions it is now possible to examine the sensitivity of 
these solutions to changes in the forcing and initial 
conditions. 

C. SENSITIVITY EXPEBIMENTS 

Using the steady state solutions from case 1 as initial 
conditions two additional experiments are carried out to 
determine the sensitivity of the planetary waves to 'errors* 
in the forcing or initial conditions. In the first experi- 
ment errors are introduced in the forcing field. In fact a 
1005? error was introduced by turning off the forcing alto- 
gether. The model is then integrated out to 20 days. Figs. 

4.17 and 4.18 contain the harmonic dials for the above inte- 
gration. Note that for the most part the once quasi- 
stationary solutions are now almost purely transient since 

the phase through the entire range of phase. In the previ- 
ously observed pattern the stationary plus transient nature 
of the waves was indicated by a change of phase and ampli- 
tude which occured in a particular quadrant of the dial. 
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The estimated phase speeds obtained from Figs 4.17 and 4.18 
show these transient waves to again be moving at the speed 
of the mean wind. The fact that the waves are pure tran- 
sient now is not surprising since the heating is no longer 
present to fix the wave to a particular location. 

In the second experiment errors are introduced into the 
rotational part of the steady state wind by reducing those 
components by 2055. The model is again integrated out to 20 
days. Figs. 4.19 and 4.20 contain the harmonic dials for 
this integration. From these figures one can see that the 
only wave which is moving eastviard is the temperature wave. 
All the other component waves are moving westward. The 
temperature wave has a phase speed of about 12 °/day at the 
equator which is about 2555 slower than the mean wind. It is 
evident that there are a number of high frequency components 
present in the other component waves which are not present 
in the temperature wave. Since there is no simple circular 
pattern for the wind and pressure waves it is difficult to 
obtain the period of the transient part of the solution. 
The exact phase speed of these wave is not important. The 
important thing to note here is that the waves are moving 
slowly westward. 
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Figure 4,17 Same as Fig 4.2 for the Case with No Heating. 
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Figure 4.18 Same as Fig 4.3 fcr the Case with No Heating. 
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Figure 4.19 Same as Pig 4.2 for the Case 
with Botational Hind Errors. 
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Figure 4.20 Same as Fig 4.3 for the Case 
with Botational Bind Errors. 
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V. DISCDSSIOH ANE CONCLOSIONS 



The results of the previous section are best discussed 
in terms of normal mode theory. A number of investigators 
(Geisler and Dickinson, 1976, Daley et al, 198 1, and 
Kasahara 1976) have used this method to examine the tran- 
sient part of planetary waves. The method basically 
consists of solving the linearized set of eguations 
describing oscillations of a stratified resting atmosphere 
on a spherical earth by separating the eguations into a 
latitudinal structure equation and vertical structure equa- 
tion. The vertical structure equation determines the 
vertical modes of oscillation. The separation 
constants (eigenvalues) from the vertical structure equation 
can be used to find the eigenfunctions of the latitudinal 
structure equation. The latitudinal structure equation for 
a given eigenvalue is identical to a fundamental equation of 
a free surface, if the eigenvalue is replaced by the uniform 
depth of the fluid. This equation was first obtained by 
Laplace in his study of free oscillations of shallow water 
over a rotating sphere. For this reason the horizontal 
structure equation is often referred to as Laplace's tidal 
equation and the eigenvalues are referred to as equivalent 
depths. The mode associated with the largest equivalent 
depth is referred to as the external or barotropic mode. 
This is because this mode has a vertical profile of hori- 
zontal divergence with the same sign throughout the atmos- 
phere. The other modes associated with smaller equivalent 
depths are referred to as internal modes. These modes are 
oscillatory in the vertical. For a given equivalent depth, 
zonal wavenumber and meridional mode there exist three sets 
of eigenfrequencies and eigenfunctions. Two are eastward- 
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and v€stward-propaga ting gravity waves and the third is a 
westward propagating Rossby-Haurwitz type wave. The phase 
speeds of these waves are related to equivalent depth, meri- 
dional mode and zonal wavenumber. For the Rossby type wave 
the fastest phase s peed occurs with the largest equivalent 
depth and smallest meridional mode and zonal wavenumber. 
The period of this gravest, symmetric zonal wavenumber one 
Rossby mode is almost exactly 5 days. There is observa- 
tional evidence for the existence of this wave (Hadden and 
Julian, 1972) and there is also evidence that some of the 
planetary wave error in numerical models is in part due to 
the spurious excitation of these large scale external Rossby 
modes (Daley et al, 1981). The following transient modes 
which were observed in the three different parts of this 
study are: 

1. the mode which moved at the speed of the mean wind 

in the case where the model atmosphere was being 
brought to a steady state, 

2. the mode which moved at the speed of the mean wind in 
the case where the model heating was turned off, 

3. the westward moving modes in the case where error was 
introduced into the rotational part of the wind. 

These modes may be explainable in terms of these external 
and internal Rossby modes. The mode which is moving at the 
speed of mean wind may be a slow westward moving internal 
mode which is overwhelmed by the mean wind. The westward 
moving modes in the last experiment are most likely internal 
modes as well but are probably associated with larger equiv- 
alent depths than those of the earlier experiments. Since 
phase speed for a given wavenumber and meridional mode 
increases with increasing equivalent depth, these waves 
would have a greater westward phase speed and would be less 
counterbalanced by the mean wind. Given a large enough 
westward phase speed one would actually see a westward 
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propagation indicated on the harmonic dials. The low 
freguency modes in the first two parts of this study were 
most likely excited by the discontinuous nature of the 
forcing. Numerical studies by Geisler and Dickinson (1976) 
indicate that even when the forcing is gradually turned on 
that low freguency Eossby modes are excited. An analytic 
study by Clark (1972) also showed that a mode which moved at 
the speed of the mean wind would be excited by a switch on 
of a vertical velocity or temperature disturbance. However, 
this mode was only one of a nunber of modes which would be 
excited and it was also found to decay with time. Thus, it 
is possible that the transient modes which were observed to 
move at the speed of the meat wind are in fact internal 
Eossby modes. The fact that internal Eossby modes might 
have been excited when errors were introduced into the rota- 
tional part of the wind is consistent with the results of 
Daley et al (1981) who found that errors in the rotational 
wind will cause a spurious excitation of these modes. 
However, it must be pointed out that since the results of 
this study were analyzed using spherical harmonics rather 
than normal modes it is impossible to confirm the above 
hypothesis. The guestion that might now be asked is why 
wasn’t the external barotropic node excited in these experi- 
ments? The studies by Geisler and Dickinson and Clark found 
that this mode would be excited. Clark found that this mode 
would in fact be the dominate mcde. The crucial difference 
between the analytic model of Clark , the numerical model of 
Geisler and Dickinson and the model in this study is that 
this model does not have a vertical wind shear while the 
aforementioned models do. This lack of vertical shear is 
the most probable reason why the barotropic mode was not 
excited. In the absence of a mean vertical shear the 
heating specified in this model will produce low pressure at 
the surface and a high pressure aloft. The vorticities 
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associated with each will be of the same magnitude so that 
the vertically averaged vorticity vanishes and there can be 
no barotropic mode response. Khen a mean vertical shear is 
present there is a differential advection affect and the 
upper and lower vorticities dc not exactly cancel, thus 
allowing the possiblity of a barotropic mode 

Given the slow phase speed of the transient modes 
excited in this study by 'errors' in the forcing and initial 
conditions it does not appear that planetary waves are very 
sensitive to errors. However, this apparent lack of sensi- 
tivity is most likely due to the simplified conditions used 
in the model rather than to the inherent nature of planetary 
waves. Specifically, the lack of wind shear prohibits the 
excitation of the fast external Rossby wave. The external 
mode has a fast phase speed ( 72 °/day) and if excited it 
could rapidly propagate throughout the region of integra- 
tion . 

The results obtained by heating in this study might also 
have been obtained by includirg topography along with a 
stronger surface friction tern. Future studies should 
investigate the sensitivity of planetary wave to errors in 
topography as well as heating. In addition, future studies 
will have to use a more realistic vertical wind profile. 
Due to the lack of a realistic wind profile, which prohibits 
the excitations of the external barotropic mode, any conclu- 
sions about the sensitivity of planetary waves to initial 
conditions and forcing must be delayed until more comprehen- 
sive studies are completed. 
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